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Abstract 

We have studied the changes in the electronic structure and the superconducting transition 
temperature T c of Mg(Bi_ x C x ) 2 alloys as a function of x with < x < 0.3. Our density- 
functional-based approach uses coherent-potential approximation to describe the effects 
of disorder, Gaspari-Gyorffy formalism to estimate the electron-phonon matrix elements 
and Allen-Dynes equation to calculate T c in these alloys. We find that the changes in the 
electronic structure of Mg{Bx_ x C x ) 2 alloys , especially near the Fermi energy Ep, come 
mainly from the outward movement of Ep with increasing x, and the effects of disorder 
in the B plane are small. In particular, our results show a sharp decline in both B and C 
Vx{y) states for 0.2 < x < 0.3. Our calculated variation in T c of Mg{B\^ x C x ) 2 alloys is 
in qualitative agreement with the experiments. 
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1 INTRODUCTION 



The changes in the electronic properties of MgB 2 upon doping with various ele- 
ments gets manifested in the changes in its superconducting behavior. An under- 
standing of these changes in the superconducting behavior of M gB 2 upon doping 
can provide useful information about the nature of interaction responsible for super- 
conductivity. Such information can then be used to prepare materials with desired 
superconducting properties. Since the discovery of superconductivity in MgB 2 , a 
substantial amount of effort has been directed towards probing the various aspects 
of superconductivity using doping. 

The effects of doping MgB 2 with various elements such as s-elements Be [1] and 
Li [2], sp-elements Al [3,4] and Si [5], and (/-elements Fe, Co, Ni [6] and others 
[7] have been studied experimentally. Most of these substitutions take place in the 
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Mg plane of MgB 2 , leaving the B plane in tact. For some dopants at particular 
concentrations, it is also found that the dopants prefer to form a separate layer [3,4] 
rather than mix with the Mg layer or the B layer. For dilute amount of dopants the 
changes in the superconducting properties, in particular the superconducting tran- 
sition temperature, T c , depends on the alloying element [1,2,3,4,5,6,7]. For larger 
concentrations, with only Mg\_ x Al x B 2 [3,4] being studied extensively, the T c is 
found to decrease due to the gradual filling of the B p x ^ -holes. 

The doping of MgB 2 with C [8,9,10,1 1,12,13] can be expected to be different from 
the ones described in the previous paragraph due to the following two reasons, (i) 
C prefers to go to the B plane instead of the Mg plane and (ii) C forms strong 
covalent-bonded layered structure similar to B layer in MgB 2 . These reasons im- 
ply that the lowering of T c in Mg(Bi_ x C x ) 2 alloys with increasing x is mainly due 
to the filling up of the p x ( y ) holes, and not due to disorder in the B plane. However, 
the x dependence of T c in Mg(B 1 _ x C x ) 2 alloys as obtained in some earlier exper- 
iments may not be accurate due to the difficulties in estimating its C content [13]. 
According to Ref. [13], the T c observed for x = 0.1 in Mg(Bi- x C x ) 2 alloys was 
22 K, which is lower than the result of Ref. [12]. 

Theoretical attempts at understanding the changes in the electronic structure and 
how these changes affect the superconducting properties of carbon-doped MgB 2 
alloys have not gone beyond the rigid-band model approach [14] with its limited 
applicability. To be able to reliably describe the effects of C doping in MgB 2 , 
we have carried out ab initio studies of Mg(B 1 ^ x C x ) 2 alloys with < x < 0.3. 
We have used Korringa-Kohn-Rostoker coherent-potential approximation [15,16] 
in the atomic- sphere approximation (KKR-ASA CPA) method for taking into ac- 
count the effects of disorder, Gaspari-Gyorffy formalism [17] for calculating the 
electron-phonon coupling constant A, and Allen-Dynes equation [18] for calculat- 
ing T c in Mg(B 1 _ x C x ) 2 alloys as a function of C concentration. We have used the 
present approach to study several other Mg_B 2 -based alloys [19,20,21]. In the fol- 
lowing, we present our results in terms of the changes in the total density of states 
(DOS), in particular the changes in the p contributions of both B and C to the total 
DOS, as a function of concentration x of C atoms. Before we describe the results 
of our calculations we provide some of the computational details. 



2 COMPUTATIONAL DETAILS 



The charge self-consistent electronic structure of Mg(B 1 _ x C x ) 2 alloys as a func- 
tion of x has been calculated using the KKR-ASA CPA method [15,16]. We parametrized 
the exchange-correlation potential as suggested by Perdew-Wang [22] within the 
generalized gradient approximation [23]. The Brillouin zone (BZ) integration was 
carried out using 1215 k— points in the irreducible part of the BZ. For DOS cal- 
culations, we added a small imaginary component of 1 mRy to the energy and 
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used 4900 k-points in the irreducible part of the BZ. The lattice constants for 
Mg(Bi_ x C x ) 2 alloys as a function of x were taken from the experimental result 
of . The Wigner- Seitz radius for Mg was larger than that of B and C. The sphere 
overlap which is crucial in ASA, was less than 10% and the maximum I used 
was l max = 2. The present calculations were carried out using the scalar-relativistic 
Schroedinger equation. The Green's function is calculated in a complex plane with 
a 20 point Gaussian quadrature. 

The electron-phonon coupling constant A was calculated using Gaspari-Gyorffy 
[17] formalism with the charge self-consistent potentials of Mg(Bi_ x C x ) 2 ob- 
tained with the KKR-ASA CPA method. Subsequently, the variation of T c as a 
function of C concentration was calculated using Allen-Dynes equation [18]. The 
parameters used in the calculation of T c using the Allen-Dynes equation, for the 
whole range of x, were /i* = 0.09, u rrns = 400 cm" 1 , and u>i n = 834 K. 



3 RESULTS AND DISCUSSION 

In this section we describe our results for Mg(B 1 _ x C x ) 2 alloys in terms of the 
changes in the densities of states, the electron-phonon coupling constant A and 
the superconducting transition temperature T c as a function of x with < x < 
0.3. In our analysis we will emphasize those aspects of changes in the electronic 
structure which directly affect the superconducting properties of Mg{B\_ x C x ) 2 
alloys, in particular its T c . Thus, we discuss the changes in the density of states as 
a function of x, in particular the p contributions of both B and C to the total DOS 
and how it affects the electron-phonon coupling constant and consequently the T c 
of Mg{B x _ x C x ) 2 alloys. 

3.1 Densities of States 

In Fig. 1, we show the total DOS of Mg(B 1 - x C x ) 2 alloys for x equal to 0.01, 0.1, 
0.2 and 0.3, calculated using the KKR-ASA CPA method. A comparison of Figs. 
l(a)-(d) shows that with increasing x the Fermi energy, E F , moves outward and the 
DOS becomes smoother. The outward movement of Ep is due to the increase in the 
number of valence electrons with the addition of C while the increase in disorder 
in the B plane of MgB 2 flattens the peaks in the DOS. In Mg{B\_ x C x ) 2 alloys the 
outward movement of E F decreases the total DOS at the Fermi energy, N(E F ), in 
particular for x ~ 0.3, N(E F ) is a minimum. Generally, a decrease in N(E F ) leads 
to a lowering of T c , it is expected that with the increase in C concentration the T c 
of Mg(Bi- x C x ) 2 alloys should decrease. 

It is known that in MgB 2 the holes in the p x ^ bands couple very strongly to the 
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Figure 1. The calculated total density of states (solid line) of Mg{B\- x C x ) 2 alloys for (a) 
x = 0.01, (b) x = 0.1, (c) x = 0.2 and (d) x = 0.3, respectively. The dashed vertical line 
indicates the Fermi energy. 




Figure 2. The calculated p x i y \ density of states for both B (solid line) and C (dot-dashed 
line) in Mg(Bi_ x C x ) 2 alloys at (a) x = 0.01, (b) x = 0.1, (c) x = 0.2 and (d) x = 0.3, 
respectively. The dashed vertical line indicates the Fermi energy. 



optical E 2g phonon mode [24,25,26,27], and hence a change in the p x ^ bands 
directly affects the superconducting properties of MgB 2 . Thus, to get a more de- 
tailed picture of how the changes in N(E F ) affect the superconducting properties 
of MgB 2 , we have to analyze the changes in the B p x i y ) contributions to N(E F ). 
In addition, since C is known to form strong covalent-bonded layered structures, it 
may be that in Mg(B 1 _ x C x ) 2 alloys the C p x ( y ) electrons/holes also couple to the 
optical phonons the way B p x ^ electrons/holes do. 
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Figure 3. The calculated p z density of states for both B (solid line) and C (dot-dashed line) 
in Mg(Bi- x C x ) 2 alloys at (a) x = 0.01, (b) x = 0.1, (c) x = 0.2 and (d) x = 0.3, 
respectively. The dashed vertical line indicates the Fermi energy. 

To see how the p contributions to the total DOS change in Mg(Bi- x C x ) 2 alloys as 
a function of x, we show in Fig. 2 the p x ( y ) DOS for both B and C in Mg(Bi- x C x ) 2 
alloys for x equal to 0.01, 0.1, 0.2 and 0.3. The p DOS of MgB 2 calculated with the 
KKR-AS A CPA method is in good agreement with our full-potential result [28] and 
other calculations [29,30], indicating the reliability of the present approach. From 
Fig. 2 it is clear that near E F both B and C p x (y) densities of states are similar for 
the whole range of x with C DOS being somewhat smaller than the B DOS. As 
the addition of C increases the electron count, the outward movement of E F as x 
increases fills up the p^bands of both B and C, signalling a sharp change in the 
superconducting properties of Mg(B 1 _ x C x ) 2 alloys in this range. 

The changes in the p z contribution to the total DOS in Mg(B 1 ^ x C x ) 2 alloys as a 
function of x is not very significant, as can be seen from Fig. 3 , where we have 
shown the p z DOS for both B and C in Mg(Bi- x C x ) 2 alloys for x equal to 0.01, 
0.1, 0.2 and 0.3. In the dilute limit, the p z DOS due to B in Mg(B 1 _ x C x ) 2 alloy 
is similar to the p z DOS due to B in MgB 2 [28,30,31]. In addition, these electrons 
play a minor role in deciding the superconducting properties of Mg(B 1 _ x C x ) 2 al- 
loys. 

Since the changes in the superconducting properties of Mg(B 1 ^ x C x ) 2 are dictated 
by the DOS at E F , in particular the p x ( y ) contributions of both B and C to the 
total DOS, we show in Fig. 4 the variations of the total density of states N(E F ) 
and the B and C p x ^ contributions to it. As can be seen from Fig. 4, there is 
a steep drop in the p x ( y ) DOS for both B and C between x = 0.2 to x = 0.3, 
indicating the filling up of the p x < y ) bands. Given that the p x i y ) bands are responsible 
for driving superconductivity in MgB 2 , we expect the superconducting properties 
of Mg(B 1 ^ x C x ) 2 alloys to change significantly between x = 0.2 to x = 0.3. 
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Figure 4. The calculated variation in the total (upper panel) and the p x ( y ) (lower panel) 
densities of states at the Fermi energy as a function of x in Mg(B\- x C x ) 2 alloys. In the 
lower panel, both B (filled circle) and C (filled triangle) p x ( y ) densities of states are shown. 

Here, we like to point out that the range of x as indicated above, in which the 
superconducting properties of Mg(B 1 - x C x ) 2 are expected to change significantly, 
must be viewed within the constraints inherent in the present approach. 

3.2 The superconducting transition temperature 

We have used the Gaspari-Gyorffy formalism to calculate the Hopfield parameter 
and then the electron-phonon coupling constant A of Mg(B 1 _ x C x ) 2 alloys. Our 
calculations show that the C atoms in Mg(B 1 _ x C x ) 2 alloys couple to the phonons 
as strongly as the B atoms. For example, at a; = 0.2 we find A (per atom) for B 
and C to be 0.31 and 0.24, respectively. It is interesting to note that as the C p x ( y ) 
bands fill up somewhat later than the B p x ^ bands, the calculated A of C becomes 
greater than that of B for x ~ 0.3. 

To see how the changes in the electronic properties of Mg(B 1 ^ x C x ) 2 alloys affect 
the superconducting transition temperature, we show in Fig. 5 the calculated T c as 
a function of x in Mg{B\^ x C x ) 2 alloys. The experimentally observed T c values as 
reported in Refs. [12,13] are also shown in Fig. 5. Note that for x = 0.1, Refs. 
[12] and [13] find T c to be equal to 34 K and 22 K, respectively. Our calculations 
show that the T c drops to ~ K from 20 K as x changes from 0.2 to 0.3, which is 
consistent with the filling up of the p x ( y ) bands as discussed above. The calculated 
trend in the variation of T c in Mg(B 1 _ x C x ) 2 alloys is in qualitative agreement with 
the experimental results as shown in Fig. 5. 

Some of the differences between the calculated (more accurately than the present 
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Figure 5. The calculated (rilled circle) variation of T c as a function of concentration x in 
Mg(B\- x Cx)i alloys. The experimental results of Ref. [12] (filled triangle) and Ref. [13] 
(star) are also shown. 

approach) and the observed T c in Mg(Bi- x C x )2 alloys may be due to an inaccu- 
rate determination of the C content in these alloys, as pointed out in Ref. [13]. 
In the present case, the calculation of T c incorporates two additional approxima- 
tions, (i) use of a constant phonon frequency over the whole range of x and (ii) use 
of Gaspari-Gyorffy formalism for estimating the electron-phonon matrix elements, 
which is known to underestimate the strength of the coupling. However, we expect 
the trend in the variation of T c in Mg(B 1 ^ x C x ) 2 alloys, as shown in Fig. 5, to be 
reliable. 



4 CONCLUSIONS 

We have studied the changes in the electronic structure of Mg(Bi^ x C x ) 2 alloys as 
a function of x with < x < 0.3 using the KKR-ASA CPA method. We find that 
the changes in the electronic structure of Mg(B 1 ^ x C x )2 alloys , especially near the 
Fermi energy Ep, come mainly from the outward movement of Ep with increasing 
x and the effects of disorder in the B plane are small. In particular, our results show 
a sharp decline in B and C p x ( y ) states for 0.2 <x< 0.3. Our calculated variation 
in T c of Mg(B 1 _ x C x ) 2 alloys is in qualitative agreement with the experiments. 
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